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fNJ ' Abstract 

> : 

The QCD sector of the system SANC is presented. QCD theoretical predictions for several 
processes of high energy interactions of fundamental particles at the one-loop precision 

■ level for up to some 3- and 4-particle processes are implemented. 

1 Introduction 

Oh 



Oh: 



The computer system SANC is aimed to carry out semi-automatic calculations at the one-loop 
precision level of realistic and pseudo-observables for various processes of elementary particle 
interactions to be investigated at the present and future colliders - Tevatron, LHC, ILC and 
others. 

We created the QCD environment of SANC and started to implement systematically NLO 
QCD processes filling the QCD branch of SANC tree in the same spirit as in Ref. PQ. 

We created a set of FORM [2] procedures for the analytic calculation of building blocks 
of QCD such as self-energies of quarks and gluons, vertices with virtual gluons and corre- 
sponding counter terms. These blocks are placed into the QCD Precomputation level of the 
system. The FORM programs are accessible via the same menu sequences as for QED or 
EW Precomputation. 

We consider here the results of implementation of the first processes of the QCD branch 
available in SANC. They are subdivided into 31egs and 41egs branches as one can see in Fig. HJ 
The 31egs branch contains b2q decays, namely: t —>■ W b, W — > ud, W — > cs, Z —>■ qq and 
H —>■ qq. Here b and q denote any weak boson and any quark, respectively, b denotes the 
b-quark. 

The 41egs branch contains 4f processes. For the latter there is a branch for Neutral 
Current (NC) processes that contains the Drell-Yan process qq — > t~ £ + , and a branch for 
Charged Current (CC) processes that contains the decay t — > b£ + ue and the Drell-Yan 
2f — > 2f process ud — > £ + vg. 
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Figure 1: QCD part of the SANC tree. 



The structure of these branches is the same as the corresponding structure in the EW 
sector of SANC. For each process there are three FORM modules: (FF) Form Factor, (HA) 
Helicity Amplitudes, and (BR) Bremsstrahlung. 

Once the three FORM codes for the calculation of (FF), (HA) and (BR) have been 
compiled and outputs transferred to the software package s2n.f, one can get the numerical 
results. The user guide for running SANC is given in Ref. pQ. 

We convolute the partonic sub-process cross section with quark density functions to get 
the cross section at the hadronic level. One must avoid double counting of the quark mass 
singularities, subtracting them from the density functions. 



SANC version vl .00 is accessible from servers at Dubna http://sanc.jinr.ru/ (159.93.74-10) 



and CERN http://pcphsanc.cern.ch/ (137.138.18042). 



2 QCD environment 

The QCD environment of SANC is a set of FORM procedures relevant for QCD. The basic 
procedure is QCDAlgebra.prc, which calculates color weights for process diagrams. It uses 
some common relations for T a matrices and structure constants f abc . We have changed 
several intrinsic SANC procedures like FeynmanRules .pre, MakeAmpSquare .pre, Trace. pre 
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by including gluon-quark and gluon-gluon vertices in the FeynmanRules .pre and giving color 
indices for quark bispinors. Quarks and leptons bispinors have different representations in 
SANC. 

For example, the Born amplitude of the Drell-Yan process with a charged current taken 
from procedure VirtCC4f QCD .pre is of the following form: 

1 q 2 
BornDYCC = -i — ^W- 
8 s-Ml 

x Vb(ii, pi, hi, ell) -y(ii, mu) 76(h) U(ii, p2, h2, cl2) 5(cll, cl2) 

x Ub(jj,p3,h3) 7 (jj,mu) 7 6(jj) V(jj,p4,h4). (1) 

Here U(ii, p2, h2, cl2) and Vb(ii, pi, hi, ell) are bispinors of the incoming up quark and down 
antiquark, ell and cl2 are their color indices. While Ub(jj, p3, h3) and V(jj,p4,h4) are 
bispinors of the outgoing lepton pair; 7(ii, mu) and 76(H) = I(ii) + 7 5(ii) are Dirac matrices. 
Here 

M 2 W = M 2 W - iM w T W} (2) 

where M w is the mass and Y w is the width of the W boson. 

Using this environment we build a set of precomputation files. The user can find it in the 
system. To compute the quark self-energy, one follows the sequence: 

QCD — > Precomputation — > Self — > Quark — > Quark Self 

in the QCD tree of SANC. Precomputed quark self-energies are used by FORM programs which 
calculate quark counter terms: 

QCD — > Precomputation — > Self — > Quark — > QuarkRenConst. 

After we have the one-loop amplitude of a given QCD process free from ultraviolet divergences, 
we can obtain a virtual radiative correction using the universal procedure MakeAmpSquare . pre 
to calculate the modulus squared of the amplitude. 

3 QCD radiative correction to b2q decays 

We started to work with simple processes of boson decays suitable to test the QCD envi- 
ronment of SANC. One-loop Feynman diagrams (see Fig. [2]) of these processes contain only a 
vertex with one virtual gluon and two quark legs, and corresponding QCD counter terms. 

Following the standard procedure of SANC, we calculate the virtual part of QCD corrections 
to these boson decays. 

The real part of QCD corrections is created by gluon emission from every quark leg (see 
Fig. El). 

Gluon bremsstrahlung implementation needs another set of procedures, specific for every 
process. Procedures BremWqqQCD.prc, BremZqqQCD.prc, BremHqqQCD.prc prepare the corre- 
sponding amplitudes for W — > qq, Z — > qq, H — > qq processes . Procedures Sof tWqqQCD .pre, 
Sof tZqqQCD .pre, Sof tHqqQCD .pre calculate the soft gluon bremsstrahlung contributions an- 
alytically. Procedures HardWqqQCD .pre, HardZqqQCD.prc, HardHqqQCD . pre prepare fully dif- 
ferential expressions of hard gluon bremsstrahlung contributions. These expressions may be 
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Figure 2: QCD vertex with two quark legs and self energy for quarks. 
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Figure 3: Gluon bremsstrahlung from two quark legs. 
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used by Monte Carlo generators (or integrators). The Hard procedures continue the analyti- 
cal calculation of hard gluon contributions by integrating over the angle between one of the 
quarks and the gluon and over the invariant mass of the final particles. 

Here we present QCD corrections to boson decays. Analytical expressions are too cum- 
bersome to be presented in this paper. One can access these results in the system. Here we 
cast the corrected decays width in the form 



■plloop 



-iBorn 



[1 + 5(m b ,m qi ,m q2 )} 



(3) 



where m&, m qi and the masses of the boson and two quarks and S is the correction. 

In regard to vector boson decays everything is straightforward. In the Table [T] we give 
one-loop numerical results for the function 5(my, m qi , m q2 ) in percent, where my is the mass 
of the vector bosons. 



process 


S{m b ,m qh m q2 ) 


% 


W -^cs 


5(m w ,m c ,m s ) 


+3.44 


Z ^bb 


5(m z , m b , m b ) 


+3.88 


Z — » uu 


5(m z ,rn u ,m u ) 


+3.41 



Table 1: Function S(mv,Tn q i,m q 2) in percent. 
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The well known formula for vector boson decay into massless quarks is: 

Mm 5(m H ,m q ,m q ) = ^-C f ^, C f = ~, (4) 

m q ^0 4 7T 3 

which gives 3.41% for all decays of vector bosons. One can see that our numbers are in a 
good agreement with the classic result. However, mass effects are significant. 

QCD radiative corrections to the Higgs boson decay into a quark pair have been considered 
earlier by one of the authors in Ref. [3]. There the one-loop QCD correction was presented 
keeping the masses of the outgoing quarks (Eq. 4.1) and neglecting these masses everywhere 
except in logarithms (Eq. 4.3). In the framework of SANC we reproduced these results. Here 
we give the function 8(mH, m qi , m q2 ) of the QCD correction in for the Higgs boson decay 
keeping the masses of the outgoing quarks in logarithms only: 

lim 5(m H ,m q ,m q ) = C f — 

77lq >() 7T 

The term with the large logarithm, 



(6) 



was obtained first and discussed in Ref. [I]. It doesn't violate the Kinoshita-Lee-Nauenberg 
theorem since the Higss-quark coupling constant is proportional to the quark mass. Moreover, 
resummation of these large logarithms in all orders of the perturbation theory is possible as 
suggested in Ref. [I]. 
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4 QCD radiative corrections to semi-leptonic top quark 
decay 

Here we discuss in detail the analytic calculation of QCD radiative corrections to the semi- 
leptonic mode of the top quark decay t — ► b£ + V£. 

First we consider the easier "cascade" calculation of the top decay: 

t — y b + W+ 

I — * e + + v t 

We use the formula for the process with production of an unstable particle (here it is W 
boson) given in the book [5] to treat the top decay partial width: 

1 t^M+u t — , { i ) 

1 w 

where IV is total width and Y w +_ >i+Ve is partial width of the W boson. The formula (j7j) is 
valid when the total width of the unstable particle (i.e. W boson) is much less then the mass 
of this particle. The order of magnitude estimate of the relative precision for this formula is 
T w /M w . 

There are not any QCD corrections to the partial width Y w +^+ V( of the W boson. So we 
take here this partial width in Born approximation. The QCD corrections are present only 
in the partial width T t _> bW + of the top quark. 
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decay 


t -> bW + 


W+ -> e+z/ e 


t — ► be + v e 
cascade 


pBorn 


1.5930 


0.22018 


0.16405 


plloop 


1.4764 


0.22018 


0.15204 


6,% 


- 7.32 





-7.32 



Table 2: Cascade approximation for QCD corrections to the top quark decay t — > be + 



4.1 QCD radiative correction to decay t — > 6VK + 

This mode of the top decay is treated in SANC in the same way as the b2q decays considered 
in the previous section. The virtual part of QCD corrections to the decay t — > bW + is similar 
to that of the decay W~ — ► du. One can prepare it by the procedure VirtTopWbQCD.prc 
and then calculate the gluon bremsstrahlung contribution using procedures specific for decay 
t -> bW + : BremTopWbQCD.prc, SoftTopWbQCD.prc, HardTopWbQCD.prc. In result we obtain: 



t->bW+ 



-iBorn 



where 5 is the QCD correction to the decay t — > feW^ 4 " . 

The authors of the Ref. [6] have studied the effects of the total width IV in cascade 
calculations of EW corrections to the top quark decay due to the photon emittion from the 
W boson. However here we deal with gluons, they are not emitted from the W boson, so the 
QCD correction to the decay t — > bW + does not depend on IV. Results of the narrow width 
cascade approximation for QCD corrections to the top quark decay t — > be + v e are given in 
Table H 

We see that the QCD correction to T t _+ bi + I/e is the same as the QCD correction to T t ^ bW +. 
This is so because according to the formula fl7|) we have 



r 



"pBorn 



[i + 5) r*T_ 



■e+u t 



w 



-pBorn 

1 t-*be+v t 



1 + 6). 



(9) 



4.2 One-loop QCD amplitude of decay t ^ b + l + + 

The one-loop QCD amplitude comes from the following gauge independent set of diagrams: 
The QCD part of the amplitude has the same structure as the electroweak one (see Ref. []]): 



A = V»J ^(Pi) [ +* 7m (1 + 75) F?l D ( s ) + il, (1 - 75) 72?° (a) 

+D, (1 + 75 ) F?™(s) + Dp (1 — 75) Jg^a) ] U t (p 2 ) 
1 



x ■ 



Ml 



C^(Ps)7m(1 + 7b) Kip*). 



(10) 



Here _D M and 4-momentum conservation read 



D^= (p 1 +p 2 ) f ,, P2=Pl+P3+P4 
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Figure 4: One-loop QCD diagrams of t quark decay: vertex and two counter terms. 



the invariant s is 



{P2-P1) 



(12) 



V t b is the element of the CKM matrix; U, U and V are the corresponding bispinors. 

Form factors obtained by FORM code FF (see Fig. [1]) are already free from ultraviolet 
divergences: 



a, 

47T 



Cf 



+ 



111 \ m l 



+ In 



mi 



^X(s,m 2 ,m 2 b 



where 



(3(m 2 ,m 2 b ,s) 
a s m t mb 

2n \/A(s, m 2 , mf) 



2s 



- 2 {m\ + m 2 b - s) C (-m 2 , -m\, -s; m t , m g , m b ) 



.m t ,m b ,s 



mi — mt 



2s 



In 



m 



mi 



m t) m b) s 



a s m t 

IT S 

n 3 m b 



In 



mi 



+ 



m 2 J (3(m 2 ,-s,ml) 



L (m 2 t ,m 2 b ,s 



mt 



In MM + 



m (3(m 2 ,s,m 2 b ) 



m t ,m b1 s 



(13) 



A(s, m 2 , ml) = s 2 + m\ + m b — 2sm 2 — 2sm\ — 2m\ 



2 2 

m 



a( 2 2 x v /A ( s ' m *' m ') fl / 2 2x yjK^mlm 2 



f3(m 2 t ,-s,ml) 



2 i 2 ' 

+ — s 
yf\(s,m$,ml) 

~ 2 2~~ ' 

— s — 



2 i 2 ' 



L(x,y,z) 



In 



1 + P(x,y,z) 
1 - P(x,y,z)' 



(14) 



The gluon infrared singularity is regularized by a fictitious gluon mass m 9 . The Passarino- 
Veltman function Co{—m 2 , —ml, —s;m t ,m g ,mb) is infrared singular. The lepton mass is 
neglected. 

The helicity amplitudes for this mode of top quark decay are the same as those given in [1] 
but we have to take QCD form factors -F^ CD (s), -F$? D (s), ^ D CD (s) and instead of 

EW ones. 
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4.3 Virtual QCD correction 

The three particle phase space element is 



ds \f\(s,m 2 ,m 2 b ) d cos $i ds du 



2tt 



8nm 2 



16vr 



128tt 3 



rn 



2 ' 



(15) 



where the invariant u is related to the angle between the lepton momentum pi in the 
R-frame {pz + Pa = 0) and the momentum p[ of the b quark, as follows: 



u 



- m 2 t +m 2 b -s + ^X(s,m 2 ,ml) costfj 

— L 

The angle i?j varies from to n and the invariant s varies in the interval 

m 2 < s < (m t — m b ) 2 . 

Using the formula 

dT = -^-\V tb \ 2 \A\ 2 d^ 3 \ 
2m t 



(16) 



17) 



(18) 



we obtain by the procedure VirtTop3f QCD. pre the virtual QCD correction to the differential 
width of the top quark decay ^ lrt [ ' — - . Here we give this analytic expression in the form: 



d 2 T Vht (s,u) 
ds du 



ds du 

d 2 T (s,u) as 
ds du 2n 

+ I v t 12 Gf 



^{s) - m t ^{s) - m b ^{s)] 



<tb 



M A w s 



a. 



167r 3 mf I s - Af2 I 27r 
-m t m b ^i CD (s) +m t u F^ D (s) +m b u ^ D (s) 



(19) 



where 



d 2 T (s,u) 
ds du 



v tb | 2 G2p 

167r 3 m 3 



Mi 



s-M 2 w 



(m 2 t -uj (u- m 2 b ) . 



(20) 



is the differential width of the top quark decay in the Born approximation. 



4.4 Gluon bremsstrahlung corrections 

The gluon bremsstrahlung amplitude is prepared by the procedure BremTop3f QCD. pre. 
The conservation of 4-momentum reads 

P2 = Pi +P3 +P4 +P5- (21) 

The energy of the emitted gluon is obtained in the rest system p 3 + p 4 + p 5 = 0: 



P5 



s — s 



(22) 





g(ps) 

Figure 5: Gluon bremsstrahlung diagrams of t quark decay. 



where s' is the invariant mass of the two final leptons, s' = — (p3 + Pa) 2 . It varies in the 
interval : 

mj < s' < s. (23) 

We see from ( |22J) , (IT7|) and f[2"3"j) that the energy of the emitted gluon could go to zero even 
when the invariant s goes to its minimal value mf . 

An auxiliary parameter uj separates the soft and hard gluon contributions. Gluon energy 
for the soft gluon bremsstrahlung lies in the limits 



< VI < 



(24) 



where u is arbitrary small. Therefore the soft bremsstrahlung amplitude is factorized by the 
Born amplitude and the kinematics of the soft bremsstrahlung is Born-like. The four-particle 
phase space element in this case is a product of the same three-particle phase space element 
d& 3 ' and the phase space element of the emitted soft gluon. 

By the procedure SoftTop3f QCD.prc the soft gluon contribution is obtained in the form: 



^ 2 r S oft(s,^) _ d 2 T (s,u) ^ as ^ ^ / s 
ds du dsdu * tc Soft ' 



(25) 



where the integration in 



«*) = -J di [ln(^)+ln(l-e 2 ) 



Pi 



2pi • P2 



+ 



(pi-n) 2 {pi ■ n){p 2 ■ ri) (p 2 -n) 



(26) 



over the angle of the gluon emission is performed in the same rest system p3 + Pa + p5 = = 
p2 — pi- Taken the gluon momentum p 5 as p^ = np® here we use the unit 4- vector n = (1, ri) 
where n has the direction of the vector p 5 . The standard method [7] by t'Hooft and Veltman 
was applied. 

The infrared divergences here are the same as the ones in the virtual gluon contribution 

^ lrt [ - — - but with the opposite sign. So, the sum of virtual and soft gluon contributions 

ds du 

does not contain any infrared divergences. We give here the expression of this sum integrated 
by IntcTop3f QCD.prc over the invariant u (i.e. over the angle 
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dT Virt (s) dTsoft(s) dT (s) a 
+ : — : °/ 



ds 



ds 

At 



ds 



2tt 
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■ L {m 2 ,m 2 b ,s) - 2 



2 2 

m f ,s,m h ) H + 1 



-In 
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1 

(3 (ml mis) 
1 



L 2 (m 2 t ,s,m 2 b ) +4 Li 2 I - 



L 2 (m 2 , -s, ml) + 4 Li 2 



L(mlm 2 b ,s) 
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+ TT^ oT £( m ?> s, m 2 ) - 4 



m£ _ _ 
2s 2s 
2fi(m 2 ,s,m 2 J 
+ /3(m 2 ,s,mf) 
' 2/?(m 2 ,-s,m 2 ) ^ 
1 + /5(m 2 ,-s,m 2 ) y 

L(m 2 t ,-s,m 2 b ) 



(5(m 2 t ,s,m 2 b ) 
^\(s,m 2 t ,m 



+ 



m t m b 



L(m 2 t ,m 2 b ,s) 



-I K 



2s yj\(s,m1,ml) 
2 G 2 F m b Cja s s (s — m 2 t — ml + 4m t mb) 
64:ir 3 m 2 ix 



(27) 



s-Ml 



L(m 2 t ,m 2 b ,s) 



The hard gluon contribution is produced by the procedure HardTop3f QCD.prc. The gluon 
energy p® for the hard gluon bremsstrahlung varies in the interval 



w < p® < max^). 
The four particle phase space element is 

ds ^\(s,m 2 t ,m 2 b 



(2? 



s — s') ds' d cos #i d cos $40^4. 



128vr 3 m 2 128tt 3 s 
For the hard gluon contribution the invariant s' varies in the interval: 

m 2 p < s' < s — 2uj \ fs. 



(29) 



(30) 

The separating parameter uj is arbitrary small. The sum of soft and hard gluon contribution 
to the decay width has no trace of it. 

The kinematics and choice of variables to be integrated over are illustrated in Fig. [6j The 
angle $i is between the b quark momentum (p[) and the emitted gluon momentum (ps), the 
angle $4 is between charged lepton momentum (pi) and emitted gluon momentum (p§). The 
z-axis is chosen here along the momentum p"s of the emitted gluon. 
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Figure 6: a) Lab. system; b) Rest system p% + Pa + Pb = 0; c) Rest system p 3 + p^ = 0. 



Using a formula analogous to (|T8l) we obtain the fully differential hard gluon contribution 
to the top quark decay. After integration over angles we get 



d 2 T Hard (s,s') dT (s) C f a s \ s - M* 



ds ds' 



x- 



ds it (s — s' 



Ml 



—— L(m 2 t , m 2 b , s) - 2 



+ I V t 



tb 



G F Cfa s 



Mi 



967r 3 mf 7r Is' — M 2 



7 , 3 , 



m\ + mf) 2 - (mf + m 2 b ) [~s + ~s') + 2s 2 + ss' + s 



J2 



L(m 2 ,ml,s) 



+ -yA(s,m 2 ,m^) s(m 2 t + m 2 b ) - (2s 2 + ss' + s 



i'2\ 



(31) 



Neglecting the mass of the charged lepton we obtain the QCD radiative correction to the 
decay width: 



5 



FfZ r bi+u e (! + <*)> 
-8.48%. 



(32) 



Our result is in agreement with the result of Ref. [8] (5 = —8.5%). 

Comparing the narrow width cascade approximation result in Table [2] with the complete 
result ( 1321 we see that the QCD radiative correction to the decay width obtained by the 
narrow width formula (Q is near but not equal to the complete result. 



5 QCD radiative corrections to Drell— Yan processes 

Here we present the results for the corrections to the charged (CC) and neutral (NC) current 
Drell- Yan processes, ud — > t^vi and qq — > respectively. All formulas below are shown 

at the partonic level. 
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At first we give expressions for cross sections in the Born approximation: 



CC/ 



V, 



G1 



Mi 



ud 



NC/ 



187r | s - Ml 

2 a/ 2 2 

n a p(s,m e} m e 



3 9 m% 
2 m « + 2? 



4 / ml \ * \ Ami , , 



(33) 



where s 



-(Pi 



V2) 



Pi and P2 are 4-momenta of the initial quarks; P(s,m 2 ,m 2 ] 



Amh 



; m# is the lepton mass. Here we denote 



V (s) 
VJs) 



(34) 



= Q\Q\ + Q q Qi lxz(s) + x * z (s)} v q v e + | xz(s) | 2 + if r 

_ - V {s)-2\ Xz {s)\"{vl + lf){ir)\ 

v q = lf ] - 2Q q sin 2 9 W , v t = lf ] - 2Qi sin 2 9 W . 

The Z/j propagator ratio Xz( s ) with s-dependent (or constant) Z width is given in [9]. 

The one-loop QCD amplitude of the charged current Drell-Yan process is similar to ffTU]) 
given in the previous section and the corresponding amplitude of the neutral current Drell- 
Yan process has another form. However, we neglect the terms proportional to the masses 
of the initial quarks and therefore the virtual QCD corrections of both processes calculated 
in corresponding procedures VirtCC4f QCD. pre and VirtNC4f QCD. pre are proportional to 
corresponding Born cross sections. 

The gluon bremsstrahlung amplitudes both for CC and NC processes are prepared by 
the procedure Brem4f QCD .pre. Here we give the sum of soft and virtual gluon contributions 
which does not contain any infrared divergences. 



j.CC 
'Virt 



cc 

Soft 



- NC | -NC 
°Virt ~r °Soft 




(35) 



(36) 

where m q , m u , m^ are quark masses. Hard and soft gluon bremsstrahlung contributions are 
calculated in procedures Sof tCC4f QCD .pre, HardCC4f QCD . pre for CC and Sof tNC4f QCD . pre, 
HardNC4f QCD. pre for NC, respectively. We present here expressions for hard bremsstrahlung 
with extracted splitting function: 

1 s 2 + s' 2 

'} 



uo Hard 



ds' 



Hard 



ds' 



2-7T S 2 S — S 1 

a„ 1 s 2 + s' 



In 


^) 


+ In 




\ m uJ 










In 


\ m2 q ) 


- 1 









- 2 



m 



(37) 
(38) 
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where s' is the invariant mass of final leptons. 

One-loop radiative corrections contain terms proportional to logarithms of the quark 

masses, In ( ] . They come from the initial state radiation contribution including vir- 
\m 2 J 

tual, soft and hard gluon emission. These terms are in agreement with the prediction of 
the renormalisation group approach, see i.e. Ref. [10J. In the case of hadron collisions these 
logarithms have been already taken into account in the parton density functions (PDF's). 
Therefore we have to apply a subtraction scheme to avoid the double counting. Linearization 
of the subtraction procedure is done as described in the Ref. [TT] . 

In order to have the possibility to impose experimental cuts and event selection procedures 
of any kind, we use a Monte Carlo integration routine based on the Vegas algorithm [T2]. In 
this case we perform a 4(6)-fold numerical integration to get the hard gluon contribution to 
the partonic (hadronic) cross sections. To get the one-loop QCD corrections we add also the 
contributions of the soft gluon emission and the virtual QCD loop. The cancellation of the 
dependence on the auxiliary parameter Q in the sum is observed numerically. 

For numerical evaluations we take the same set of input parameters as the one given in 
Ref. [13] . Below we consider the obtained distributions of corrections 8 for Drell-Yan CC and 
NC processes, where 

da l-loop 

° = —r~5 1- 

The total QCD corrections to the CC and NC cross-sections are small, however the shape of 
distributions changes greatly. 
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Figure 7: Transverse momentum P-j and transverse mass Mj of neutrino-lepton pair distri- 
bution for Drell-Yan CC. 

Fig. EH) and Fig. [7b) present the charged lepton transverse momentum pj and neutrino- 
lepton pair transverse mass Mj distributions of QCD correction 5 for Drell-Yan CC process, 
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where Mj = y2p^pj (1 — cos and </^„ is the angle in the transverse plane. The be- 
haviour of distributions at left side is an edge effect - near the cat we have a deficit of the 
hard gluon emission. Below the resonance the correction is small. Near to the resonance pj 
distribution becomes negativ and after the resonance rapidly grows up-to 250 per sent. But 
the My distribution remains small in an wide range. 

At Fig. [Hk) we show the lepton pair transverse momentum distribution of QCD correction 
S for Drell-Yan NC process and at Fig. [Hb) - lepton pair invariant mass distribution. Their 
behaviour is similar to the CC case. 
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a) b) 
Figure 8: Transverse momentum P t and invariant mass Mu of lepton pair distributions for 
Drell-Yan NC. 



Comparison of QCD and EW distributions was discussed in reports on workshop |14j . 
We are going to develop a Monte Carlo event generator to describe the Drell-Yan processes 
in realistic conditions. 



6 Summary 

In this paper we described the first steps of creation of the QCD branch in the SANC system. 
We developed environment for calculations of QCD processes. Then we tested successfully 
this environment by simple calculations of b2q decays. We implemented corrections to top 
quark decays and observed that complete one-loop result differs from the narrow width cas- 
cade approximation result. Drell-Yan NC and CC processes were implemented into SANC, 
and we created the Monte Carlo integrators to study them. We are going to develop uni- 
fied event generator for Drell-Yan processes, which would include QCD and EW corrections 
simultaneously. 
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